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In  this  thesis,  we  investigated  the  performanee  of  wireless  sensors  in  caseaded  multi-hop 
relay  network  in  Underwater  Aeoustie  Communieation  (UAC)  to  see  the  optimum  data  rate  that 
ean  be  transmitted  through  the  ehannel  with  a  minimum  bit  error  rate.  We  explored  the 
experimental  evaluation  and  simulation  result  for  the  mutual  information  between  transmitted 
signal  and  reeeived  signal  in  aquatie  ehannel.  We  optimized  the  mutual  information,  finding 
optimum  channel  capacity  for  a  cascaded  multi-hop  channel  network.  Our  results  revealed  that 
the  placement  of  multiple  sensor  nodes  over  a  long  range  in  underwater  communication 
improves  the  channel  capacity  in  a  significant  way.  A  relay  system  network  with  multiple  sensor 
nodes  at  different  distances  reduced  the  propagation  loss,  information  loss,  multipath  fading  etc. 
and  minimizes  the  bit  error  rate  (BER). 

Since  UWC  is  bandwidth  limited,  we  studied  the  relationship  of  capacity  with 
transmission  power,  bandwidth  and  carrier  frequency.  In  comparison  with  the  free  space  wireless 
communication,  the  underwater  acoustic  communication  suffered  from  the  limits  of  the  less 
available  bandwidth,  multipath  effect,  and  the  complex  noise  caused  by  the  underwater  acoustic 
channel.  Therefore,  we  explored  a  relay-aided  UAC  in  order  to  enhance  reliability  in  long 
distance  underwater  communication  and  defy  the  challenging  problems  such  as  frequency 
dependent  signal  attenuation,  long  propagation  delay,  and  information  loss.  The  result  illustrated 
better  performance  in  terms  of  improvement  in  channel  capacity,  minimum  BER,  and  the 
increase  in  mutual  information  compared  to  traditional  direct  link  UAC.  In  addition,  effects  of 
various  system  parameters  on  channel  capacity  were  considered  such  as  source  to  destination 
distance,  transmit  power  allocation,  sensor  node  location,  carrier  frequency  etc.  To  demonstrate 
the  merits  of  above  discussed  network  model,  we  simulated  its  BER  performance.  Comparison 
was  made  with  traditional  direct  link  systems  to  reveal  the  communication  reliability  of  new 
network  model. 

We  focused  on  channel  capacity  as  a  function  of  SNR  for  different  link  distances.  Our 
contribution  on  maximization  of  mutual  information  between  transmitted  information  and 
received  information  was  done  by  determining  the  optimal  carrier  frequency  and  node  distance. 
Thus,  optimal  frequency  was  derived  in  terms  of  transmission  distance  and  power  allocation.  Our 
findings  proved  that  the  optimal  frequency  decreases  with  an  increase  in  range.  Similarly, 


capacity  of  the  relay  ehannel  inereased  with  an  inerease  in  the  number  of  hops  for  varying  link 
distanees. 


Chapter  1 


Introduction 

In  this  thesis,  we  investigate  the  performance  of  acoustic  sensor  networks  in  underwater 
wireless  communication  in  terms  of  mutual  information,  channel  capacity,  bit  error  rate,  and 
information  loss  based  on  the  transmission  power,  distance,  carrier  frequency  and  the  bandwidth 
of  the  channel.  A  placement  of  multiple  acoustic  sensors  in  a  multi-hop  relay  form  is  needed  in 
numerous  applications  where  data  transmission  has  to  be  accomplished  beyond  short  distances. 
This  network  is  effective  since  the  bandwidth  of  an  underwater  acoustic  communication  is 
severely  limited  and  decays  with  increase  in  distances.  It  is  advantageous  to  accomplish  such 
transmission  using  sensors  in  a  multi-hop  relay  form  keeping  constraints  such  as  transmission 
rate,  transmission  delay,  Signal-to-Interference  and  Noise  Ratio  (SINR)  under  consideration.  In 
particular,  we  consider  a  communication  scenario  where  a  certain  number  of  bits  have  to  be 
transmitted  over  a  distance  d.  Our  result  analyzes  the  optimal  number  of  sensor  nodes,  N  to  use 
over  a  link  distance  d  to  increase  the  mutual  information  between  transmitted  data  and  received 
data  in  terms  of  link  capacity,  delay,  information  loss,  and  bit  error  rate. 

This  work  is  supported  by  the  grant  titled,  “Information  -  Driven  Doppler  Shift  Estimation 
and  Compensation  Methods  for  Underwater  Wireless  Sensor  Networks,”  (Department  of 
Defense  Research  Grant  W911NE-1 1-1-0144)  which  considers  fundamental  research  problems 
[11].  The  proposal  centralizes  innovative  information  driven  approaches  to  make  a  decision  for 
sensor  collaboration  related  to  the  network  information  constraints. 

Different  goals  of  the  above  mentioned  project  can  be  summarized  as  [1 1]: 

1)  To  analyze  and  improve  the  non-data  aided  methods  for  Doppler  shift  estimation  and 
compensation  (squaring  time  phase  recovery  method,  power  spectrum  method,  and  partial 
autocorrelation  method); 

2)  To  determine  the  necessary  quantity  and  deployment  strategy  of  sensors  in  a  given 
region,  and  to  provide  a  required  security  level; 

3)  To  minimize  the  probability  of  false  alarm  and  the  bit  error  rate; 
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4)  To  improve  decision-making  on  target  detection  with  sensor  collaboration  in  the 
context  of  blind  Doppler  shift  estimation  and  detection,  and 

5)  To  improve  the  ability  of  a  distributed  sensor  system  to  detect  intruders  (targets)  and 
determine  how  the  network  wireless  sensor  design  is  affected. 

Our  thesis  concentrates  on  improving  mutual  information  between  two  data  signals  in 
terms  of  minimum  bit  error  rate,  maximum  data  rate  communication,  and  channel  performance. 
We  focus  on  sensor  placement  in  multi-hop  to  minimize  the  bandwidth  consumption  and 
mitigate  the  risk  of  acoustic  link  failure  and  information  loss.  In  order  to  reduce  the 
communication  load  for  each  sensor  we  develop  new  distributed  algorithms  for  underwater 
communication  similar  to  those  developed  in  wireless  communication  in  [14].  The  advantage  of 
this  sensor  placement  will  help  in  fulfilling  the  increasing  demand  for  reliable  maximum  data 
rate  wireless  communication  links  to  accommodate  the  wide  range  of  underwater  application. 
We  also  define  algorithms  that  reduce  the  significant  challenges  to  the  development  of 
underwater  wireless  systems  such  as  complex  channel  noise,  bandwidth  limitations,  and 
frequency  and  distance  dependent  signal  attenuation  [5]. 

According  to  Shannon  (1948),  the  mathematical  theory  of  channel  capacity  is  based  on  the 
notion  of  maximizing  the  mutual  information  between  the  input  and  the  output  of  a  system.  The 
growing  demand  for  underwater  wireless  communication  gives  great  importance  to  determining 
capacity  limits  of  acoustic  channels  for  different  network  models.  These  capacity  limits  show  the 
maximum  data  rates  that  can  be  transmitted  over  wireless  channel  with  small  error  probability. 
Shannon  defines  channel  capacity  as  the  channel’s  mutual  information  maximized  over  the  all 
the  possible  input  distributions.  UWC  impose  many  constraints  that  affect  the  design  of  wireless 
network  such  as  path  loss,  transmission  distance,  and  energy  absorption  by  water,  long 
propagation  delays,  channel  dispersion,  Doppler  shift,  and  interference.  However,  after  the  wide 
development  in  wireless  digital  communications  for  the  underwater  environment,  it  is  important 
to  improve  the  mutual  information,  data  rate,  channel  capacity,  information  loss,  probability  of 
error  etc.  of  existing  systems.  In  our  thesis,  we  modeled  our  sensors  as  a  multi-hop  relay  network 
to  maximize  mutual  information  for  point  to  point  Gaussian  channels  between  transmitter  and 
receiver  in  the  underwater  environment. 
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Our  contribution  is  on  UWSN  design  for  long  range  aeoustic  communication  placing  an 
optimum  number  of  sensor  nodes  at  an  optimum  distanee  for  certain  average  transmission  power. 
A  simple  mathematical  model  assessing  the  dependence  of  an  aeoustic  communication  channel 
capaeity  on  the  distance  is  given  in  [3],  where  the  aeoustie  path  loss  experienced  by  a  signal  of 
frequency / propagating  in  a  horizontal  distance  d  was  modeled  as  A(d,  where  k  is 

the  spreading  faetor  and  a(f)  the  absorption  eoefficient.  The  eapaeity  ealeulation  is 
aecompanied  with  two  types  of  bandwidth,  i.e.  a  3  dB  design  prineiple  and  an  optimal  eapacity- 
maximizing  prineiple. 

Similarly,  maximizing  the  mutual  information  with  an  upper  bounded  average  power  yields 
the  same  ehannel  capacity  [6].  The  algorithm  in  [6]  focuses  on  relationship  between  power, 
transmission  band,  distance  and  eapaeity  for  the  Gaussian  noise  scenario  in  a  closed  form.  This 
model  estimates  for  power  consumption,  band  edge  frequency  and  bandwidth  as  functions  of 
distance  and  capacity  required  for  a  data  link.  This  approximate  model  stems  from  an 
information  theoretic  analysis  that  takes  into  aceount  a  physical  model  of  aeoustic  propagation 
loss  and  both  physical  and  ambient  noise.  It  shows  that  the  power  and  band  of  operation  can  be 
adjusted  to  maximize  mutual  information  to  a  certain  capacity  level. 

Both  carrier  frequency  and  bandwidth  selection  are  important  in  UWSN  design  for  multi-hop 
relay  networks.  For  long  range  eommunication,  frequeney  has  to  be  lowered  and  the  available 
bandwidth  to  work  will  be  just  few  kHz.  We  focus  on  non-coherent  channels  for  maximizing  the 
mutual  information  between  transmitter  and  reeeiver.  We  study  the  transmission  distanee,  data 
rate,  and  channel  noise  to  get  the  maximum  of  the  digital  information  embedded  in  the  received 
signal  to  maximize  the  mutual  information. 

In  this  thesis,  we  applied  a  multi-hop  relay  sensor  node  by  dividing  the  total  link  into 
multiple  hops,  where  a  relay  aeoustic  node  is  employed  at  each  hop.  The  relay  node  reeeives  the 
signal,  regenerates  it,  and  passes  it  on  to  the  next  hop,  until  the  final  destination  is  reaehed.  The 
question  of  interest  to  our  design  is  how  exaetly  does  this  approach  gives  improvement  in  terms 
of  maximizing  the  mutual  information,  total  power  consumed,  energy  per  bit,  overall  eost  and 
delay.  We  use  an  information  theoretic  tool  to  accomplish  the  assumption  that  the  eapaeity  of 
relay  ehannel  is  equal  to  the  capacity  of  each  of  its  hops,  and  are  shown  to  inerease  with  the 
number  of  relays.  Here,  we  show  that  relaying  also  helps  to  reduce  the  total  transmission  power 
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and  its  benefits  are  even  more  pronouneed  in  view  of  energy  per  bit  savings.  We  consider  system 
optimization  in  the  light  of  minimizing  both  bit  error  rate  and  the  information  loss. 

Multi-hop  relay  acoustic  sensor  network  channels  are  useful  in  system  design  in  which 
transmission  rate  is  optimized  according  to  the  number  of  hops.  The  result  makes  use  of  a  basic 
acoustic  path  loss  model,  thus  providing  guideline  for  the  design  of  a  general  relay  acoustic  link. 
We  calculate  the  optimal  carrier  frequency  and  input  signal  PSD  to  maximize  mutual 
information  and  improve  channel  capacity. 

1.1.  Thesis  Structure 

To  understand  the  underwater  acoustic  communication,  its  channel  capacity,  information 
loss,  and  signal  power  allocation;  one  requires  knowledge  of  information  theoretic  tools,  digital 
communication,  and  signal  processing.  Also,  we  have  extensively  used  the  Matlab  programming 
in  system  design  and  simulation  of  sensors  placement  to  maximize  the  data  rate,  mutual 
information,  and  channel  capacity.  It  is  also  very  important  to  explain  a  few  details  in 
appreciable  the  thesis. 

The  rest  of  the  thesis  is  structured  is  structured  as  follows: 

Chapter  2  introduces  our  sensors  network  model  for  sensors  placement  in  multi-hop  relay  form. 
Chapter  3  provides  the  detail  about  the  optimal  frequency  and  distance  for  underwater 
communication  along  with  the  bit  error  rate  and  delay  that  incurs  due  to  multi-hop  relay  network. 
Chapter  4  explains  the  proposed  model  of  multi-hop  relay  network  in  underwater 
communication  valid  from  Information  theoretic  point  of  view  too.  It  also  briefly  discusses 
entropy,  joint  Entropy,  information  loss,  and  mutual  Information. 

Chapter  5  discusses  about  the  modulation  techniques  for  the  scenario  of  multi-hop  relay 
network  in  underwater  wireless  communication. 

Chapter  6  does  the  simulation  and  analytical  study  for  optimal  distance,  frequency,  transmitted 
power  as  a  function  of  number  of  hops  in  an  Underwater  Acoustic  Channel. 

Chapter  7  contains  the  main  observation  and  results  from  experiment  performed  at  swimming 
pool.  We  studied  the  mutual  information  and  information  loss  from  the  observation  of  bit  error 
rate  due  to  channel  noise,  distance,  and  carrier  frequency  in  underwater  acoustic  communication. 
Chapter  8  provides  the  conclusion  and  scope  for  future  work. 
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1.2  Thesis  Significance,  Contribution  and  Findings 

This  main  scope  of  thesis  is  to  maximize  the  mutual  information  between  transmitter  and 
receiver  in  long  range  communieation  in  an  underwater  channel  for  high  data  rate  transmission 
using  a  multi-hop  relay  network.  It  foeuses  on  the  mathematieal  algorithm  for  long  range 
horizontal  aeoustie  eommunieation  for  optimal  frequency,  power  alloeation,  distanee  and 
available  bandwidth.  We  use  direet  sequenee  spread  speetrum  (DSSS)  and  PN-sequenee 
modulation  as  better  methods  of  modulation  that  is  robust  to  noise  and  interference. 

This  researeh  eontributes  to  Dr.  Cotae’s  UDC  researeh  titled  “Information  -  Driven  Doppler 
Shift  Estimation  and  Compensation  Methods  for  Underwater  Wireless  Sensor  Networks” 
(Researeh  Grant  W911NE-1 1-1-0144),  whieh  is  to  analyze  and  develop  non-eoherent 
oommunieation  methods  at  the  physieal  layer  for  target  traeking  and  use  the  information  theory 
tools  to  prediet  the  next  target  position  by  proeessing  the  data  information  eolleeted  from  a 
eollaborative  wireless  distributed  sensor  network. 
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Chapter  2 


Underwater  Acoustic  Communication  Channel 

In  this  chapter,  we  provide  a  baekground  review  on  eoncepts  used  throughout  the  thesis.  We 
review  path-loss,  effeet  of  noise  and  interferenee,  propagation  delay,  and  long  range  aeoustie 
eommunieations.  We  review  information  theoretie  tool  terminologies  such  as  mutual 
information,  entropy,  joint  entropy,  eonditional  entropy,  and  ehannel  eapaeity. 

2.1  Underwater  Channel  Model 

An  underwater  eommunieation  ehannel  is  very  eomplex  and  eommunieates  through 
propagation  of  sound,  ealled  aeoustie  eommunieation.  Underwater  aeoustie  sensor  networks  are 
then  formed  by  aeoustieally  conneeted  bottom  sensor  nodes,  autonomous  underwater  vehieles, 
and  surfaee  stations  that  serve  as  gateways  and  provide  radio  eommunieations  link  to  on-shore 
stations.  Underwater  aeoustie  sensor  networks  eonsist  of  sensor  nodes  and  vehicles  deployed 
underwater  and  networked  via  the  aeoustie  links  to  perform  eollaborative  monitoring  tasks. 
However,  aeoustie  channels  impose  many  eonstraints  that  affeet  the  design  of  UWC  systems. 
These  are  charaeterized  by  path  loss  that  depends  on  both  the  transmission  distanee  and  the 
signal  frequeney.  The  signal  frequeney  determines  the  absorption  loss,  whieh  inereases  the 
distanee  as  well,  eventually  imposing  a  limit  on  the  available  bandwidth. 


Table  1.  Underwater  Acoustic  Communication  System  Ranges  [3] 


Distance 

Range(Km) 

Bandwidth  (KHz) 

Short 

0.1 

>100 

Very  short 

0.1-1 

20-60 

Medium 

1-15 

10 

Long 

15-100 

2-6 

Very  long 

1000 

<1 
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Table  I  shows  typical  bandwidth  of  the  underwater  channel  for  different  ranges.  The  table  is 
generated  using  3dB  bandwidth  module  for  relation  between  relative  SNR  and  frequency  from 
Figure  11.  This  shows  that  long-range  systems  that  operate  over  kilometers  may  have  a 
bandwidth  of  only  a  few  kHz,  while  a  short-range  system  operating  over  few  meters  may  have 
more  than  a  hundred  kHz  bandwidth.  This  is  true  because  carrier  frequency  is  higher  for  short 
range  acoustic  communication  which  will  give  a  wide  bandwidth  to  transmit  the  signal. 
Likewise,  for  long  range  acoustic  communication,  carrier  frequency  is  lowered  and  available 
bandwidth  is  small. 

2.1.1  Acoustic  Transmission  Loss 

Acoustic  path  loss  depends  on  both  signal  frequency  and  distance.  The  dependence  is  a 
consequence  of  absorption,  i.e.  acoustic  energy  is  transferred  into  heat  energy  during  sound 
propagation  from  one  point  to  another. 

A(d,/)  =  d“a(/)^  (2.1.1) 

where  d  =  total  transmission  distance  (km) 

/  =  carrier  frequency  (kHz) 

a(/)  =  absorption  coefficient  (dB/km) 

a  =  spreading  factor 

In  addition,  signal  experiences  a  spreading  loss  which  also  increases  with  distance.  Spreading 
loss  refers  to  the  energy  distributed  over  an  increasingly  larger  area  due  to  the  regular  weakening 
of  a  sound  signal  as  it  spreads  outwards  from  the  source.  The  overall  transmission  loss  that 
occurs  in  UW  channel  over  a  transmission  distance  at  a  signal  frequency /is  given  by  [3]; 

TL  =  a  10  log  d  +  dlO  loga(/)  (2.1.2) 

where  a  is  spreading  factor  (a  =2  for  spherical  spreading,  a  =1  for  cylindrical  spreading,  and  a 
=  1.5  for  the  practical  spreading),  d  is  transmission  distance  in  meters.  In  general,  for  shallow 
water  channels,  cylindrical  spreading  is  assumed  (oc  =\)  while  for  deep  water  channels  spherical 
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channels  is  assumed  (a  =2).  And  10  log  a(f)  is  the  absorption  coefficient  expressed  using  the 
Thorp’s  formula,  which  gives  a(f)  in  dB/km  for  f  in  kHz  as  follows  [3]: 


101oga(/)  =  0.11 


1+f 


■  +  44 


4100  +/ 


+  (2.75)10-^^  +0.003 


(2.1.3) 


Figure  1  Absorption  coefficient  vs.  frequency 

2.1.2  Noise 

Noise  in  an  underwater  channel  affects  communication  system  by  altering  signal  strength 
and  reducing  channel  capacity.  The  noise  in  acoustic  channel  constitutes  of  four  different  sources: 
turbulence,  shipping,  waves,  and  thermal  noise.  Turbulence  noise  influences  on  very  low 
frequency  region,  /  <  10  Hz.  Noise  caused  by  distant  shipping  is  dominant  in  the  frequency 
region  10  Hz  -100  Hz,  and  it  is  modeled  through  the  shipping  activity  factor  s,  whose  value 
ranges  between  0  and  1  (for  low  and  high  activity).  Surface  motion,  caused  by  wind-driven  waves 
{w  is  the  wind  speed  in  m/s)  is  the  major  factor  contributing  to  the  noise  in  the  frequency  region 
10  Hz  -100  kHz  (which  is  the  operating  region  used  by  the  majority  of  acoustic  systems).  Finally, 
thermal  noise  becomes  dominant  for  /  >  100  kHz  [2]. 
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The  noise  in  an  aeoustie  underwater  ehannel  under  discussion  is  complex  additive  colored 
(frequency  dependent)  Gaussian  noise  with  zero  mean  and  power  spectral  density  (PSD)  N^if) 
that  decays  linearly  on  the  logarithmic  scale  in  the  frequency  region  1-20  kHz.  We  use  the  same 
approximation  as  given  in  [2]  for  the  PSD  of  noise  in  an  acoustic  channel 

101ogiV^(/)  =  iVi-77log(/)  (2.1.3) 

for  the  positive  constant  iVi  =  50  dB  re  p  pa  [24]  and  // =  18  dB/decade  [21].  Obviously,  the 
noise  PSD  of  a  narrowband  underwater  network  depends  on  the  carrier  frequency  but  it  is  flat 
over  the  communication  band.  Hence,  the  noise  effect  in  underwater  wireless  communication  is 
carrier  frequency  dependent. 

2.1.3  Channel  Capacity  and  Signal-to-Noise  Ratio 

The  channel  capacity  strictly  depends  in  transmission  distance,  carrier  frequency,  and 
bandwidth  in  an  underwater  acoustic  communication.  The  channel  is  time  variant  and  the 
ambient  noise  is  assumed  to  be  Gaussian.  Channel  capacity  is  the  information  rate  over  the 
acoustic  link  which  is  bounded  by  the  Shannon  capacity  defined  as: 

C  =  (BW)  log2  (1  +  SNR)  (2. 1 .4) 

where  BW  is  bandwidth  of  the  system  in  kHz  and  SNR  indicates  the  signal  to  noise  ratio. 

C  =  BWlog2(l  +  ^)  (2.1.5) 

where  P  =  Power  of  signal  (Watts  or  Volts  ) 

N  =  Power  of  noise  (Watts  or  Volts  ) 

It  can  be  seen  from  equation  2.1.5  that,  for  a  source  with  constant  power  and  limited 
bandwidth,  the  capacity  of  a  channel  is  highly  dependent  on  the  noise  power.  Since,  noise  and 
attenuation  is  a  function  of  distance  and  frequency,  using  A(d,f)  and  PSD  iVo(/),  the  signal-to- 
noise  ratio  (SNR)  at  receiver,  at  distance  d  and  frequency  / for  a  transmitted  power  of  P  is  given 
by  [6]: 


SNR 


P 

BWMdjjNjn 


(2.1.6) 
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Combining  equation  2.1.4  and  2.1.6,  we  get  the  channel  capacity  for  an  acoustic  link  for  optimal 
link  distance  and  frequency  as  below  [6]; 

C  =  (BW)  log,(l  +  SNR)  =  iBW)  log,  (l  +  (2.1.7) 

2.2  Network  model:  Multi-hop  relay  sensor  network 

The  multi-hop  relay  networks  have  been  extensively  studied  in  terrestrial  environments  to 
provide  reliable  communications  with  extended  transmission  distance.  However,  due  to  peculiar 
characteristics  of  the  UAC  channel,  the  benefits  of  multiple  sensor  nodes  in  UAC  needs  to  be  re¬ 
evaluated,  and  requires  special  design  and  modulation.  In  our  research,  a  multi-hop  relay  setup  is 
considered;  where  all  relays  locate  on  line  and  are  equally  spaced  with  same  transmit  power.  We 
investigate  the  maximum  data  rate  transmission  and  mutual  information  for  a  relay  node  through 
analytical  and  numerical  results.  Accordingly,  we  compare  performance  of  channel  capacity  with 
the  traditional  direct  link  system,  and  affecting  factors  of  the  system  capacity  like  power 
allocation,  link  distance,  available  bandwidth,  information  loss,  and  bit  error  rate  (BER). 

First  we  discover  an  end  to  end  direct  link  systems  BER,  and  information  loss  corresponding 
to  an  optimum  carrier  frequency  and  signal  bandwidth.  Eater,  same  observation  is  performed  for 
multi-hop  relay  network  and  mutual  information  gain  in  UAC  is  evaluated.  The  affecting  factors 
of  channel  performance  are  studied,  which  includes  the  data  rate,  distance  from  source  to 
destination,  and  carrier  frequency.  We  optimize  the  system  for  performance  evaluation  with 
respect  to  maximum  mutual  information  and  channel  capacity  for  desired  low  BER  and 
minimum  information  loss. 

We  consider  a  relay  link  as  shown  in  Figure  1.  The  total  link  distance  d  is  divided  into  N 
multi-hops,  such  that  the  length  of  each  hop  is  d/N.  The  link  is  designed  so  that  the  signal-to- 
noise  ratio  at  the  input  to  each  relay  equals  the  target  value  SNRq.  This  design  ensures  a  fair 
comparison  between  a  multi  hop  and  a  single-hop  link.  In  what  follows,  we  shall  use  the  capacity 
in  the  optimal  sense,  i.e.  the  bandwidth  and  the  power  is  defined  as  those  that  maximize  the 
channel  capacity  under  the  constraint  of  fixed  total  power  needed  to  ensure  the  desired  SNR. 
With  each  type  of  bandwidth  selection,  there  is  an  associated  transmission  power  required  to 
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achieve  pre-speeified  SNR  at  the  reeeiver.  In  partieular,  it  is  shown  in  [6]  that  for  a  given 
required  SNR,  the  following  relationships  hold; 


5(x)  =  Boid/do)-^  =  bd-^ 

C(x)  =  Co{d/do)~^  —  cd~^ 

P(x)  =  P^{d/doY  =  pd'P  (2.2.1) 


Eaeh  of  these  quantities  is  thus  eharaeterized  by  two  model  parameters:  the  sealing  faetor  h  for 
bandwidth,  c  for  capacity,  p  for  power  and  exponent  y,  (p  respectively. 


surface 


source 


relay 


relay 


destinauon 


Dotlom 


Figure  2  Illustration  of  Multihop  network. 

2.3  Conclusion 

In  this  ehapter,  underwater  aeoustie  communieation  channel  model  and  parameters  are 
reviewed.  We  discuss  the  dependeney  of  link  distance  and  carrier  frequeney  in  attenuation  and 
ambient  noise.  Similarly,  our  network  model  and  approaehed  aeoustie  sensors  network  topology 
and  parameters  relating  to  UWSN  are  explained. 
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Chapter  3 


Communications  Channel  and  Algorithm 

In  this  chapter,  we  mathematieally  define  and  derive  the  algorithm  for  multi-hop  relay 
eommunieations  ehannels  where  we  analyze  an  optimum  plaeement  of  sensor  in  long  range  of 
oommunieation  link.  We  studied  an  algorithm  that  will  optimize  the  BER  and  information  loss 
for  given  link  distanee  and  proposed  relay  network. 

3.1  System  Optimization 

Using  an  expressions  (2.2.1)  for  N  multi-hops  and  transmission  link  distanee  d,  the  ehannel 
eapaeity  of  eaeh  hop  is  determined  as  C(d/N)  in  terms  of  bits  per  seeond,  and  the  eapaeity  of  eaeh 
node  link  is  equal  to  [7]: 


C„(d)  =  C(d/iV)  (3.1.1) 

A  similar  relationship  ean  be  established  for  the  bandwidth, 

S„(d)  =  S(d/iV)  (3.1.2) 

where  S(d/iV)  is  the  bandwidth  of  eaeh  hop  measured  in  kHz. 

The  bandwidth  effieieney,  defined  as  follows  a  reverse  trend  whieh  means  it 

deereases  as  the  hop  distanee  beeomes  shorter.  However,  the  absolute  value  of  the  bandwidth 
effieieney  does  not  ehange  mueh  for  the  range  of  values  d  and  number  of  hops  N. 


single  hop 


iy- 


o— 


multiple  hops 
O - O - O 


TX  I- 
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Figure  3  Non-cooperative  linear  multi-hop  relay 

Similarly,  the  power  needed  to  transmit  over  a  one  hop  is  equal  to  P(d/iV)  measured  in  Watts 
or  Volts  ,  and  the  power  needed  to  span  the  entire  link  is: 


P„(d)  =  NP(d/N) 


(3.1.3) 
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(3.1.4) 


=  2nH  X  lQ^  ‘^^°S(l/N)+^^^^^+0.0002(l/N)+OASNRo+6m5 

SNRq  —  SNR  at  transmitter  node 
H  =  Water  depth 


p(0  =  SNR  (/,/)^(/,/)iV(/)(5M/)  (3.1.5) 

As  it  ean  be  expeeted,  less  power  is  required  with  deployment  of  numerous  nodes  or  using 
sensor  relays.  While  sueh  situation  is  typical  of  any  wireless  channel,  the  heavy  dependence  of  the 
capacity  on  the  number  of  relays  is  particular  to  the  UAC.  This  fact  bears  influence  on  the  energy 
per  bit  needed  to  transmit  over  the  distance  d.  If  the  transmission  is  accomplished  over  A  hops,  the 
energy  per  bit  is: 


_  P„(d)  _  NP^jd) 

Cpj(^d)  Cpj(^d) 


(3.1.6) 


The  equation  above  shows  that  how  the  energy  per  bit  Epj(d)  decreases  with  number  of  hops. 
The  benefits  of  multi-hopping  are  most  evident  in  this  metric,  as  it  combines  the  power  reduction 
and  the  capacity  (bandwidth)  increase  simultaneously  as  the  given  link  is  divided  into  more  hops. 


After  this,  optimum  operating  frequency  is  selected  from  algorithm  such  that  there  is  minimum 
path  loss  and  BER  is  minimal.  The  signal-to-noise  ratio,  SNR(d,f)  which  is  equivalent  to 
maximizing  the  channel  capacity,  leads  to  the  maximization  of  Saif).  From  chapter  2,  equation 
2.1.7, 


SNR^d.f) 

where  SdCf)  is  the  power  spectral  density  of  received  signal  in  Watts  per  Hertz. 

Given  the  required  rate  for  the  link  R,  the  bandwidth  BW  and  the  frequency  /  we  obtain  A(d,  /) 
as: 
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A(d.n  = - -  (3.1.8) 

(BW)No(f)(2  /bw  -  1) 

3.2.  Bit  Error  rate  Simulation 

Bit  error  rate  performance  of  a  receiver  is  a  figure  of  merit  that  allows  different  designs  to 
be  compared  in  a  fair  manner.  The  first  step  to  analyze  the  mutual  information  and  capacity  of 
system  is  to  find  BER.  Using  the  bit-error  probability  and  Shannon’s  communication  theory,  the 
upper  bound  of  mutual  information  is  obtained  both  analytically  and  in  simulations. 

Bit  error  rate  testing  with  Matlab  is  very  simple,  but  does  require  some  prerequisite 
knowledge.  BER  testing  requires  a  transmitter,  a  receiver,  and  a  channel.  We  begin  by  generating 
a  long  sequence  of  random  bits,  provide  as  input  to  the  transmitter.  The  transmitter  modulates 
these  bits  onto  some  form  of  digital  signaling,  which  we  send  through  a  simulated  channel.  We 
simulate  the  channel  by  adding  a  controlled  amount  of  noise  to  the  transmitted  signal.  This  noisy 
signal  then  becomes  the  input  to  the  receiver.  The  receiver  demodulates  the  signal,  producing  a 
sequence  of  recovered  bits.  Einally,  we  compare  the  received  bits  to  the  transmitted  bits,  and 
tally  up  the  errors.  BER  performance  is  usually  depicted  on  a  two  dimensional  graph.  The 
ordinate  is  the  normalized  signal-to-noise  ratio  (SNR)  expressed  as  Eb  /No:  the  energy-per-bit 
divided  by  the  one-sided  power  spectral  density  of  the  noise,  expressed  in  decibels  (dB).  The 
abscissa  is  the  BER,  a  dimensionless  quantity,  usually  expressed  in  powers  of  ten.  The 
achievable  data  throughputs,  and  the  reliability  of  an  underwater  acoustic  communication 
system,  as  measured  by  the  BER,  must  be  determined  to  suit  the  bandwidth  limitations  and 
distortions  of  underwater  acoustic  channels.  The  BER  tolerance  of  about  10“^  makes  it  a  viable 
model  for  poor  quality  band-limited  underwater  channels. 

Due  to  filtering  and  other  delay-inducing  operations,  for  typical  of  most  receivers,  there  is  an 
offset  between  the  received  bits  and  the  transmitted  bits.  Before  we  can  compare  the  two  bit 
sequences  to  check  for  errors,  we  must  first  determine  this  off-set.  One  way  to  do  this  is  by 
correlating  the  two  sequences,  then  searching  for  the  correlation  peak.  Suppose  our  transmitted 
bits  are  stored  in  vector  Tx,  and  our  received  bits  are  stored  in  vector  Rx.  The  received  vector 
should  contain  more  bits  than  the  transmitted  vector,  since  the  receiver  will  produce  meaningless 
outputs  while  the  filters  are  filling  and  flushing.  If  the  length  of  the  transmitted  bit  vector  is  Ltx 
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and  the  length  of  the  reeeived  bit  veetor  is  Lrx,  the  range  of  possible  offset  is  between  zero  and 
Lkx  ~  Ltx  ~1.  We  ean  find  the  offset  by  performing  a  partial  eross-eorrelation  between  the  two 
vectors. 


Figure  4.  BPSK  Simulation  for  BER  v/s  SNR  of  AWGN 

The  resulting  vector  is  a  partial  cross-correlation  of  the  transmitted  and  received  bits,  over 
the  possible  range  of  lags:  0:  Lrx  -  Ltx  ~1.  We  find  the  location  of  maximum  value  of  correlation 
which  tells  us  the  offset  between  the  bit  vectors.  Since  Matlab  numbers  array  elements  as  1 : 
instead  of  as  0:  A^-1,  we  need  to  subtract  one  from  the  index  of  the  correlation  peak. 

Once  we  know  the  offset  between  the  transmitted  and  received  bit  vectors,  we  are  ready  to 
calculate  the  bit  errors.  For  bit  values  of  0  and  1,  a  simple  difference  reveals  bit  errors.  Wherever 
there  is  a  bit  error,  the  difference  between  the  bits  will  be  +  1,  and  wherever  there  is  not  a  bit 
error,  the  difference  will  be  zero.  Using  Matlab,  we  calculate  the  error  vector  from  the 
transmitted  bit  vector,  Tx,  and  the  received  bit  vector,  Rx,  along  with  the  offset. 

The  error  vector  contains  non-zero  elements  in  the  locations  where  there  were  bit  errors.  We 
need  to  tally  the  number  of  non-zero  elements,  since  this  is  the  total  number  of  bit  errors  in  our 
simulation.  Therefore,  we  calculate  the  total  number  of  bit  errors,  Te,  from  the  error  vector. 
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Each  time  we  run  a  bit-error-rate  simulation;  we  transmit  and  reeeive  a  fixed  number  of  bits. 
We  determine  how  many  of  the  reeeived  bits  are  in  error,  and  then  eompute  the  bit-error-rate  as 
the  number  of  bit  errors  divided  by  the  total  number  of  bits  in  the  transmitted  signal. 

3.3  Delay  in  multi-hop  network 

Delay  in  UWC  is  defined  in  terms  of  transmission  delay  and  propagation  delay.  Delay  is 
viewed  as  yet  another  form  of  eost  assoeiated  with  multi-hop  sensor  node  network.  The  total 
transmission  delay  is  important  metrie  in  the  design  of  sensors  network  model  for  a  long 
transmission  distance.  In  digital  communication  system,  it  is  beneficial  that  relays  regenerate  a 
packet,  rather  than  simply  amplifying  and  forwarding  it.  This  way,  noise  aeeumulation  along  the 
link  is  prevented,  and  the  performance,  as  measured  by  the  overall  BER  stays  elose  to  that  of 
single  hop  link  with  the  same  SNR.  In  particular,  the  BER  of  an  N  -hop  link  is  given  by  the  first 
order  approximation  [7]; 

Pe,N(SNRo)  =  1  -  (1  -  Pe,^(SNRo))^  «  NPe,^(SNRo)  (3.3.1) 

where  SNRg  is  the  input  signal-to-noise  ratio  to  eaeh  relay,  and  Pg  ^  is  the  BER  of  eaeh  hop. 

The  approximation  normally  holds,  since  the  link  is  designed  such  that  Pg j^^SNRg)  «  1. 
Hence,  the  overall  BER  Pg  j^  is  only  scaled,  whieh  is  a  minor  degradation  for  a  practical  number 
of  hops.  However,  the  proeess  of  regeneration  requires  that  the  entire  paeket  be  reeeived  before  it 
is  passed  on  the  next  node,  whieh  introduees  a  delay.  We  foeus  for  simplieity  on  the  basie  delay 
that  oeeurs  due  to  regeneration,  not  taking  into  aeeount  the  time  spent  in  repeated  transmissions. 


Ti 


*■  ■ 

delay  per  hop 


Figure  5.  Delay  between  Relays 
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In  Figure  4,  the  delay  incur  in  regenerating  a  packet  at  the  n*  relay  node  is  illustrated.  If  a 
packet  contains  L  number  of  bits  (binary  modulation  can  be  assumed  without  loss  of  generality), 
then  its  duration  is: 

=  L/R  (3.3.2) 

where  Ti  is  duration  in  seconds  and  R  is  a  bit  rate  measured  in  bits  per  second.. 

Assume  a  total  link  distance  d  is  divided  equally  into  N  nodes.  Then  the  propagation  delay 
between  each  node  will  be: 


Pa  =  d/NCs 


(3.3.3) 


where  Pa  =  propagation  delay  between  each  node  (seconds) 

Cg  =  speed  of  sound  in  an  underwater  (normally  1500  m/s). 


Packet  regeneration  introduces  an  additional  delay  equal  to  the  packet  duration  T^.  Hence, 
delay  per  hop  is  the  packet  duration  plus  propagation  delay,  and  the  total  delay  that  packet 
experiences  over  A  hop  is  [7]: 


=  N 


-  d  L- 


d  L 
—  +  -N 
c.  R 


(3.3.4) 


Therefore,  the  overall  delay  increases  with  N,  and  this  fact  imposes  a  limit  on  number  of 
relays.  We  can  assess  the  effect  of  number  of  relays  into  two  cases.  First,  transmitter  rate  can  be 
fixed  regardless  of  the  number  of  hops  where  delay  will  be  linearly  proportional  to  the  number  of 
nodes.  Second,  the  transmission  rate  is  allowed  to  vary  with  the  node  distance,  so  that  it  can  be 
optimized  according  to  the  link  capacity.  In  this  case,  the  total  delay  is  no  longer  a  simple  linear 
function  of  number  of  nodes.  Instead,  it  depends  on  the  number  of  nodes  through  transmission 
rate  R. 
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Figure  6  End-to-end  delay  of  multi-hops  relay  network 


In  a  practical  system,  transmission  rate  is  proportional  to  available  link  bandwidth,  whieh 
depends  on  the  distance  between  nodes.  Ideally,  from  information  theoretic  point  of  view,  error 
free  transmission  can  be  accomplished  at  the  rate  equal  to  channel  capacity.  In  either  ease, 
available  transmission  rate  becomes  a  function  of  number  of  hops  N.  By  increasing  the  number 
of  hops,  more  packets  contribute  to  overall  delay  by  their  duration  However,  each  packet 
duration  becomes  shorter  that  accomplishes  transmission  faster  over  a  shorter  hop. 

3.4  Conclusion 

In  this  chapter,  some  insight  for  system  optimization  is  diseussed.  Optimum  distance  and 
frequency,  bit  error  rate  and  delay  for  transmission  of  large  number  of  bits  is  also  presented  and 
discussed. 
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Chapter  4 


Information  Theoretic  Analysis  in  Maximizing  Channel  Capacity 

In  this  chapter,  we  go  through  maximizing  mutual  information  and  ehannel  eapaeity 
analysis  from  information  theoretie  point  of  view.  The  basie  performanee  measure  is  ehannel 
eapaeity  that  determines  the  maximum  data  rate  that  ean  be  supported  with  an  arbitrarily  small 
error  probability. 

4.1  Mutual  Information  and  Channel  capacity 

Aeeording  to  Shannon  (1948),  ehannel  eapaeity  is  defined  as  ehannel’s  mutual  information 
maximized  over  all  the  possible  input  distribution.  Channel  eapaeity  is  an  upper  bound  of  the 
mutual  information  that  gives  the  maximum  data  transmission  rate  that  ean  be  reliably 
transmitted  through  a  eommunieations  ehannel.  The  information  ehannel  eapaeity  of  a  diserete 
memoryless  ehannel  is  [20]; 

C  =  maXpf^x^  I{X]  Y)  (4.1.1) 

This  quantity  is  defined  as  the  maximum  of  the  average  mutual  information  I{X;  Y)  between 
input  information  X  and  its  eorresponding  output  information  Y  of  the  ehannel  over  a  ehoiee  of 
distribution  ofX.  Mutual  information  satisfies  1{X)  T  )  =  l(Y  ;  A)  >  0. 

Considering  X  as  the  input  and  Y  as  the  diserete  output/deteeted  output  of  the  system,  the 
mutual  information  is  given  by  [9]; 

1{X)  Y)  =  H(Y)  -  H(Y\X)  (4.1.2) 

Entropy  is  also  ealled  self-information  beeause  /  (A;  A)  =  H(^X). 

In  other  way,  the  mutual  information  I(X;  Y)  is  the  relative  entropy  between  the  joint 
distribution  and  the  produet  distribution  [20]. 

;a;y')  =  ^^p(x,y)log^g^  (4.1.5) 

X  y 
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Figure  7.  Graphical  representation  of  Conditional  Entropy  and  Mutual  Information  ofX  and  Y 


4.2  Markov  Chain  and  Data  Processing  Inequality 

We  use  the  data  processing  inequality  to  show  that  information  loss  can  be  reduced 
forming  the  sensor  nodes  as  a  Markov  chain.  For  data  transmitted  through  X  ^  Y  ^  Z,  the 
conditional  distribution  of  Z  depends  only  on  Y  and  is  conditionally  independent  of  X.  Specifically 
the  joint  probability  function  of  forX,  Y  and  Z  can  be  written  as: 

p(x,y,z)  =  p(x)p(y|x)p(z|y)  (4.2.1) 

We  use  above  mentioned  theory  to  maximize  mutual  information  between  random  variable  X 
and  Z  by  deploying  another  sensor  node  Y  in  between  the  propagating  distance  of  X  and  Z.  This 
way,  sensor  node  Y  becomes  a  hop  in-between  transmission  distance  X  andZ.  Doing  so  reduces 
transmission  distance,  increases  bandwidth,  and  data  rate.  In  similar  way,  creating  N  nodes 
channels  to  transmit  data  from  one  point  to  another  affected  BER,  information  loss,  mutual 
information,  and  channel  capacity  of  the  system. 

Each  sensor  nodes  in  the  network  is  made  a  primary  transmitter  with  its  own  probability 
distribution  and  bit-error  probability  p.  The  objective  of  our  thesis  is  to  organize  the  sensors  in  a 
relay  network  system  to  optimize  the  transmitting  signal  power  along  long  distance  in  an 
underwater  communication. 
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Figure  8.  Cascaded  Channel  (Multi-hop)  sensor  network  model 


In  this  sensors  network  model,  we  re-transmit  the  received  signal  to  another  sensor  node 
exactly  as  they  are  received.  Thus,  in  order  to  obtain  the  input-output  statistics  of  the  cascade,  we 
need  only  to  consider  the  signal  as  one  sample  at  a  time.  Let  the  information  X  of  the  first 
transmitter  has  its  corresponding  probability  p(x).  The  information  X  will  transmit  down  the  first 
channel  and  will  be  received  as  7/  by  the  first  receiver  sensor  node,  Y2  by  the  second  receiver 
sensor  node,  and  finally  as  7  by  the  last  receiver  sensor  node.  During  the  observation  of 
downstream  transmission  of  information  data  from  X  to  7,  each  channel  will  have  different 
conditional  probability  distribution  since  the  channel  state  is  changing  with  time  i.e.  for  the 
channel  p(yt|yt-i)  gives  the  probability  distribution  of  the  samples  7. 

The  channel  capacity  for  discrete  time  Gaussian  channel  with  signal  power  ‘P’  for  bandwidth 
limited  performance  is  given  by  [21], 


(4.2.3) 


However,  by  cascading  and  using  data  processing  inequality  [20]  we  have,  for  7  ^  7  ^  1/  7, 


7(7;  7)  </(!/;  7) 

We  can  find  the  max  of  both  sides, 

C  =  maxp^x)  7(7;  V)  <  C"  =  maXp^^)  7(77;  V)  =  C(p^x)) 

Finding  capacity  of  the  product  channel  of  the  Markov  chain  [4],  we  obtain: 

7i  ^  7i  ^  72  ^  72 

7(7i,72;7i,72)  =  77(7i,  72)  -  77(7,  72|7i,72) 

=  77(7i,72)  -77(7i|7i,72)  -77(72|7i,72) 

=  77(7i,  72)  -77(7  |7i)  -77(72  172) 

<  77(7i)  +  77(7)  -  77(7 |7i)  "  77(7172) 
7(7i,72;  7,  7)  =  7(7i;  7)  +  KXp,  7) 


(4.2.5) 


(4.2.6) 
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Therefore, 


C  =  maXp(^i,,2)K^i.^2;  Yi'Yi) 

<  maXp(^ip,2)K^i;  Yi)  -  "^at^p(xi,x2)K^2.>2) 

=  maXp^xi)KXi,  >i)  +  maXp^x2)KX2-,  >2) 

C  <  C]^  +  C2 

4.3  Information  Loss 

Information  loss  (L)  between  two  different  sensors  network  models  is  ealeulated  from 
eonditional  entropy  of  input  and  output  data  [5]. 

L{X  ^Y)=  H{X\Y) 

H{X\y)  =  lyPyiy)  [-IxPx\Yix\y)\og(P^^y(ix\y))]  (4.3.1) 

Similarly, 

/f(r|x)-^p(Xi)/f(r|x  =  x) 

XiEX 

=  ^  p(Xj)  H(p)  =  H(p)  (4.3.2) 

XiEX 

Similarly,  to  caleulate  the  information  loss  in  between  the  easeaded  ehannel  system,  we 
individually  estimate  the  conditional  entropy  for  each  channel  communication. 

Intuitively,  in  cascaded  channel,  output  Y  of  input  data  X  is  directly  used  as  input  to  transmit 
to  Z  .This  gave  us: 

L(Z  ^  Z)  >  L(Z  ^  T)  +  L(Y  ^  Z)  (4.3.2) 

4.4  Mutual  Information  over  Cross  correlation 

Mutual  information  is  an  ideal  metric  for  our  sensor  placement  problem  in  multi-hop 
because  it  results  in  minimizing  the  conditional  entropy  of  our  object  state.  In  long  range 
underwater  communication,  placement  of  multiple  sensors  in  hop  increases  the  mutual 
information  among  transmitter  and  receiver  reducing  bit  error  rate,  information  loss,  and  large 
propagation  delay. 
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Mutual  information  identifies  and  quantitatively  eharacterizes  relationship  between  data  sets 
that  are  not  deteeted  by  eommonly  used  linear  measures  of  correlation.  Given  two  data  sets  as 

input  and  output  of  channel,  {X}={x\,  x2,  . xN}  and  {7}  ={yl,  y2, . yN},  their  mutual 

information,  I{X,  Y),  is  the  average  number  of  bits  of  {X}  that  can  be  predicted  by  measuring  {7} . 
Mutual  information  relationship  is  symmetrical,  I{X,  Y)  =  I{Y,  X).  Analyzing  observational  data, 
calculation  of  the  mutual  information  occurs  in  two  contexts: 

i)  Identification  of  nonlinear  correlation, 

ii)  In  the  investigation  of  causal  relationships  with  directed  mutual  information. 

Calculation  of  mutual  information  is  important  in  the  characterization  of  causal  relationships 
between  two  information  data  set.  Mutual  information  captures  complex  association  between 
variables  that  are  missed  by  standard  correlation  measures.  By  definition,  a  correlation  measures, 
either  linear  or  nonlinear,  quantifies  the  degree  of  correlation  between  {X}  and  {7}  under  their 
respective  definitions,  but  it  does  not  necessarily  identify  causal  relationships  in  the  sense  of 
identifying  which  variable  drives  the  other,  if  indeed  such  a  relationship  exists.  However,  the 
principle  of  the  mutual  information  maximization  follows  naturally  from  the  expected 
uncertainty  minimization  criterion  in  a  Bayesian  filtering  framework  [22]. 

4.5  Conclusion 

In  this  chapter,  information  theoretic  tools  like  entropy,  channel  capacity,  mutual 
information,  joint  entropy,  and  conditional  entropy  are  discussed.  We  relate  our  sensor  network 
mode  of  multi-hop  with  respect  to  information  theory  knowledge  and  also  mention  why  we 
choose  mutual  information  over  correlation  to  study  the  performance  of  network  model. 
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Chapter  5 


Modulation  Techniques 

In  this  chapter,  we  discuss  some  of  the  modulation  teehniques  that  will  give  a  better  performance 
and  optimum  information  reeeiving  teehniques  for  underwater  aeoustie  eommunieations.  In  this 
thesis,  we  study  two  different  types  of  modulation  and  use  them  in  data  eommunieation  proeess. 
PPM  is  the  first  one.  This  modulation  teehnique  is  used  by  the  vendor  (Desert  Star  Systems).  The 
basie  idea  of  PPM  is  that  the  position  of  eaeh  pulse  in  relation  to  the  position  of  a  reeurrent 
referenee  pulse  is  varied  by  each  instantaneous  sampled  value  of  the  modulating  wave.  In  our 
experiment,  we  use  m-sequenee  modulation  and  Binary  Phase  Shift  Keying  (BPSK)  modulation, 
so  we  foeus  on  these  two. 

5.1  PN-sequence  modulation 

Maximal  length  sequenees  are  generated  by  using  linear  feedbaek  shift  register.  In  [12]  it 
shows  that  shift  registers  result  a  binary  output  that  is  periodie  with  period,  iV  =  2”  —  1,  where  n 
is  the  degree  of  a  polynomial.  For  n  -stage  shift  register,  the  total  number  of  non-zero  states  will 
be  2”  —  1  all  having  period  of,  iV  =  2”  —  1  .Shift  registers  sequence  having  maximum  possible 
period  for  a  n  -stage  are  ealled  maximal-length  sequenees  or  m-  sequenees.  Maximum  length 
sequenees  are  very  useful  in  spread  speetrum  technology  because  of  following  properties: 

1.  A  maximal  length  sequenee  eontains  one  more  one  than  zero.  The  number  of  ones  in  the 
sequenee  is  -  (iV  -f  1). 

2.  The  modulo-2  sum  of  an  m-sequenees  and  any  phase  shift  of  the  same  sequenees  is 
another  phase  of  the  same  m-sequenee  (shift-and  -add  property) 

3.  If  a  window  of  width  r  is  slid  along  the  sequence  for  N-shifts  eaeh  r-tuple  exeept  the  all 
zero  r-tuple  will  appear  exaetly  onee. 

The  periodie  autoeorrelation  function  is  mathematically  expressed  as: 
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(5.1.1) 


f  1  for  k  —  IN 
fork  ^  IN 

where,  I  is  any  integer  and  N  is  the  sequenee  period. 

In  our  thesis,  we  generate  the  PN-sequenee  for  the  eode  length  of  1024  by  ehoosing  the  10**' 
order  polynomial. 

5.2  Direct  Sequence  Spread  Spectrum  (DSSS)/  Binary  Phase  Shift  Keying 
(BPSK) 

The  DSSS  proeess  is  performed  by  multiplying  a  earrier  frequeney  and  a  pseudo-noise  (PN) 
digital  signal.  The  PN  eode  is  modulated  onto  an  information  signal  using  modulation  teehniques 
binary  phase-shift  keying  (BPSK).  For  a  DSSS  system,  the  spreading  rate  is  inereased  with 
higher  number  of  ehips  per  data  symbol,  whieh  will  lead  to  an  improvement  of  the  SNR  per 
symbol  in  additive  white  Gaussian  noise  (AWGN).  For  time -varying  ehannels,  the  gain  depends 
on  the  stability  of  the  ehannels  environment  that  wave  propagates  in.  The  gain  of  the  reeeiver  is 
redueed  if  the  ehannel  ehanges  eonsiderably  during  one  symbol  period.  This  results  in  a  net  loss 
when  the  length  of  the  spreading  eode  is  inereased  past  a  eertain  point. 

Comparison  of  DSSS  and  FHSS  based  on  performanee  depends  on  the  eharaeteristies  of  an 
aeoustie  ehannel.  Aeeording  to  simulation  done  in  [10],  the  primary  limitation  is  the  time 
variability.  The  DSSS  performs  very  well  when  eoupled  with  a  ehip  rate  equalizer.  The  variation 
in  the  ehannel  gives  the  limitation  on  the  performanee  of  the  reeeiver  at  high  spreading  rates.  The 
FHSS  is  more  vulnerable  to  the  Doppler  shift  effeet  beeause  the  signals  are  transmitted  in  narrow 
bands,  but  more  robust  to  multiple  aeeess  interferenee  (MAI)  than  DSSS.  FHSS  has  a  higher  bit 
error  rate  than  DSSS,  though  FHSS  has  simple  reeeivers  and  gives  robustness  to  the  near-far 
problem  espeeially  when  used  with  eonvolution  eoding  and  soft-deeision  Viterbi  deeoding,  thus 
simplifying  the  power  eontrol  funetionality.  However,  performanee  for  FHSS  is  limited  in 
frequeney-seleetive  ehannels. 

In  DSSS  the  reeeived  signal  is  multiplied  with  the  original  pseudo  noise  eode  and  integrated 
over  the  period  on  the  data  symbol,  also  known  as  de-spreading.  The  de-spreading  will  de- 
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correlate  the  multiple  aeeess  interferenee  made  by  multiple  signals.  The  problem  with  this 
method  is  that  it  has  a  very  high  eomplexity  and  ean  be  a  problem  to  implement  in  small  sensor 
with  limited  power  [3]. 

Even  though  there  are  advantageous  and  disadvantageous  for  both  direct  sequence  spread 
speetrum  (DSSS)  and  frequeney  hopping  spread  spectrum  (FHSS)  the  DSSS  has  the  best 
performance  in  underwater  acoustic  communication.  Therefore  DSSS  is  ehosen  as  the  spread 
speetrum  teehnique  for  the  simulation  in  this  study. 


cos(2Tr  + 


Figure  9.  BPSK  Transmitting  channel  model 


cU-f,) 


Output 

data 


Figure  10.  BPSK  Receiving  channel  model 

The  advantages  of  direet  sequenee  spread  speetrum  (DSSS)  mentioned  above  are  good 
arguments  for  using  DSSS  in  underwater  acoustic  communication.  In  addition,  its  performance 
in  simulations  is  also  quite  good  and  with  some  optimization  the  results  ean  be  better.  The 
ehallenge  in  DSSS  is  to  have  good  time  synehronization  between  the  transmitted  eode  spreader 
and  the  reeeived  eode  de-spreader.  To  reduee  the  problem  of  noise,  better  fdters  ean  be  used  in 
the  reeeiver.  The  disadvantage  with  using  spreading  eodes  is  that  the  bandwidth  is  inereased  and 
this  can  be  a  problem  in  band-limited  situations.  To  deal  with  multiple  users  and  noise,  the  DSSS 
is  used  with  Gold  codes  as  spreading  codes.  These  spreading  eodes  work  for  detecting  the  correet 
symbols. 
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Chapter  6 


Simulation  and  Analysis  Study 

In  this  chapter,  we  present  an  analytieal  and  simulated  result  of  systems  ehannel  eapaeity  and 
funetion  of  proposed  multi-hop  sensor  network  in  long  range  communieation.  We  study  the 
mathematieal  algorithmic  result  of  transmission  model,  eapaeity  bound  and  optimum  frequeney 
for  link  distanee.  Sinee  we  aim  to  maximize  the  throughput  eapaeity  of  our  sensor  network  by 
the  plaeement  of  multi-hop  sensor  nodes,  in  the  following,  we  deseribe  the  system  eapaeity, 
transmission  data  rate,  and  an  optimal  distanee. 


Let  X  is  the  total  number  of  bits  of  a  paeket.  Let  Ti  represent  the  paeket  transmission  time  of 
a  sensor  multi-hop.  By  the  law  of  large  numbers,  the  throughput  eapaeity  of  the  network  is  [8]: 


C  = 


lim 

n^oo 


nX 

yn  ^ 
^t=l  ^  I 


(6.1.1) 


We  assume  the  length  of  data  paeket,  X,  is  fixed;  therefore  in  order  to  maximize  the 
throughput  eapaeity  of  the  network,  E[Ti]  needs  to  be  minimized.  Thus,  the  design  objeetive  of 
our  system  is  to  minimize  the  expeeted  time  that  a  node  needs  to  transmit  a  paeket  to  following 
node  whieh  we  eall  it  the  expected  paeket  transmission  time.  The  paeket  transmission  time  is 
inversely  proportional  to  the  eommunieation  rate  between  transmitter  and  reeeiver.  Sinee  an 
underwater  acoustie  ehannel  is  time  variant  and  different  random  network  parameters  such  as  the 
distanee  between  the  transmitter  and  reeeiver  and  ehannel  noise  affect  the  reliable 
eommunieation  rate,  the  paeket  transmission  time  is  random.  So  the  transmission  model  is  based 
on  the  Shannon  eapaeity  bound  for  a  point  to  point  Gaussian  channel. 

We  assume  that  the  noise  variance.  No,  is  the  same  anywhere  in  the  network.  Let  the  signal 
power  reeeived  by  a  reeeiver  be  P  and  the  available  bandwidth  be  BW.  The  ehannel  eapaeity  is: 

C  =  BW\og2{l +  ^)  bps  (6.1.2) 


If  the  transmitter  referenee  power  is  P,  while  the  distanee  between  the  transmitter  and 
reeeiver  is  /,  by  using  (2),  the  bit  rate  between  this  transmitter  and  reeeiver  is: 
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(6.1.3) 


C  =  BW\og2(l  +  j^j  bps 

where  a  is  the  positive  eonstant  representing  path  loss. 

Similarly,  if  the  transmitter  has  M  bits  to  transmit  to  the  reeeiver,  the  data  transmission  time 


M 

T(l,P,X)= - p — seconds  (6.1.4) 


Figure  1 1 .  Target  SNR  v/s  Capacity  maximization:  bandwidth  and  capacity  scaling  factor  bo  and  Co 

As  it  ean  be  seen  in  Figure  9,  the  targeted  SNRo  for  eommunieation  between  each  single  hop 
has  direct  dependence  on  scaling  factor  of  bandwidth  and  capacity  of  the  channel.  From 
equation  2.2.1,  the  SNR  for  total  link  is  characterized  by  parameter  bandwidth,  capacity  and 
number  of  hops.  For  multi-hop  relay  setup  of  sensors  in  underwater,  the  overall  link  capacity 
and  SNR  affects  in  optimum  selection  of  bandwidth  and  capacity  calculation  of  each  hop. 


Figure  12.  Target  SNR  under  the  Capacity  maximizing  definition:  power  scaling  factor  Pj, 
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Figure  10  illustrates  the  transmission  power  parameter  as  a  funetion  of  SNR.  The  results  indieate 
that  power  sealing  faetor  and  power  exponent  exhibit  some  SNR  dependenee  that  inereases  the 
signal  power  and  increases  the  performance  of  acoustic  channel. 


Figure  13.  Relative  SNR  in  an  aeoustie  ehannel  v/s  frequeney  as  a  funetion  of  link  distanee 
Figure  11  illustrates  that  for  each  transmission  distance  /,  there  clearly  exists  an  optimal 
frequency  /  for  which  the  maximal  narrow-band  SNR  is  obtained.  Using  the  attenuation  A(l,  f) 
and  the  noise  PSD,  once  can  evaluate  the  SNR  observed  over  a  distance  /.  In  our  simulation 
result,  we  define  3  dB  bandwidth  at  the  various  range  of  communication  for  SNR(l,f). 


Figure  14.  Optimal  frequency/!/)  with  respect  to  link  distance  d 
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Figure  12  represents  the  optimal  frequeney  plotted  as  a  function  of  transmission  distance.  In 
practice,  one  may  choose  some  transmission  bandwidth  around  f(l)  and  adjust  the  transmission 
power  so  as  to  achieve  desired  SNR  level. 


Figure  15.  Capacity  of  Relay  link  v/s  number  of  hops,  N  for  varying  link  distance  d 
In  figure  13,  the  capacity  of  multi-hop  relay  link  is  demonstrated  with  respect  to  number  of  hops. 
The  capacity  CN(d)  versus  the  number  of  hops  N  is  shown  for  varying  link  distance  d.  We 
observe  that  the  capacity  increases  with  the  number  of  hops  for  a  given  link  distance.  As 
expected  the  capacity  is  lower  for  a  greater  link  distance.  This  tells,  number  of  hops  should  be  in 
consideration  for  given  link  distance  to  find  upper  bound  in  mutual  information  between 
transmitted  data  and  received  data. 


Figure  16.  Bandwidth  of  Relay  channel  v/s  number  of  hops,  N  for  varying  link  distance 
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Figure  14  shows  the  similar  trend  like  in  capaeity  versus  number  of  hops.  The  bandwidth 
availability  for  acoustic  channel  becomes  more  with  decrease  in  communicating  distance.  The 
deployment  of  multiple  sensors  to  form  multi-hop  relay  network  for  given  link  distance  gives  an 
opportunity  to  use  wide  bandwidth  for  maximum  data  rate  transfer. 


Figure  17.  Total  Transmission  Power  Vs.  the  number  of  hops,  N 
Figure  15  represents  the  relation  between  transmit  power  and  number  of  hops  for  given 
link  distance.  As  it  can  be  expected,  less  power  is  required  for  a  greater  number  of  relays. 
Therefore,  in  underwater  acoustic  channel,  number  of  relays  shows  heavy  dependence  on  the 
capacity  of  channel. 
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Chapter  7 


Experimental  Set  Up  and  Results 


After  a  preliminary  researeh  and  signal  analysis,  we  eoneentrate  our  study  to  the  main 
objective  of  our  thesis;  an  experimental  approach  for  placement  of  sensors  in  underwater  for 
minimum  BER  and  maximum  mutual  information.  Our  experimental  model  consisted  of  an 
indoor  swimming  pool,  two  pairs  of  SAM-1  acoustic  transducers  provided  by  Desert  Star 
Systems,  two  portable  computers  both  equipped  with  MATLAB  software,  signal  processing 
toolbox,  communication  toolbox,  a  hydrophone  to  measure  sound  underwater,  and  an  acoustic 
speaker  to  generate  noise.  The  data  is  transferred  using  transducers  and  is  processed  in  Matlab. 


In  our  experiment,  the  transmitted  signal  has  a  carrier  frequency  of  17  kHz  and  the  sampling 
rate  of  8000.  From  the  Figure  17,  we  can  see  that  the  power  of  noise  is  high  for  frequency  range 
of  0  to  5  kHz.  As  the  frequency  increases  we  can  observe  that  the  power  of  noise  is  less  than 
38dB.  So  if  we  are  able  to  transmit  the  signal  in  the  frequency  range  where  the  noise  power  is 
relatively  less,  we  can  obtain  better  performance.  Therefore,  we  choose  17  kHz  of  carrier 
frequency  for  our  experiments. 


Figure  18.  Acoustic  Sensors  by  Vendor  and  Experimental  set  up  in  Underwater  Communications 
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Power  Spectral  Density 


Figure  19.  PSD  of  Noise  (No(f))  experimentally  measured  in  underwater 
Figure  17  shows  that  power  of  noise  is  high  for  frequency  range  of  0  to  5  kHz.  As  the  frequency 
increases  we  observe  that  the  power  of  noise  decreases  up  to  -12  dB.  Therefore,  we  choose  to 
transmit  the  signal  in  the  frequency  range  where  the  noise  power  is  relatively  less,  and  gives 
better  performance  in  terms  of  BER  and  mutual  information  bound. 


Figure  20.  Matlab  simulation  Set  Up  for  multi-hop  network  sensor  model 
Figure  18  gives  the  insight  of  how  the  multi-hop  relay  network  is  setup.  In  Matlab  simulation,  we 
connect  an  output  from  first  hop  to  second  hop  as  an  input  and  transmit  data  with  some  delay  and 
noise.  Before  transmitting  a  data  to  second  hop,  BER  and  mutual  information  is  calculated  and 
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Figure  21  Transmitted  and  Received  Modulated  m-Sequence  Signals  (17  kFtz) 
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data  are  regenerated  such  that  noise  and  error  does  not  get  accumulated  throughout  the 
transmission. 


Table  2.  Modem  Specifics  and  Operating  Function 


Type  of  command  /specification 

Operating  function/condition  of  sensors 

### 

Sets  to  a  command  mode  for  configuration. 

We  used  S5  for  transmission  @  13bits  /sec 

Set  acoustic  transmit  data  speed 

We  used  R5  for  receiving  (^13bits/sec 

Set  acoustic  receive  data  speed 

Serial  buffer 

32  byte 

Operating  temperature 

0-70  "C 

Filter  type 

4“^  order  band  pass  filter 

Operating  frequency  for  single  channel  receiver 

33.8khz 

Operating  frequency  of  modem 

33khz-42khz  Omnidirectional 

Sonar  Data  format 

16-  bit  data  +  4  bit  checksum  per  word 

Modulation 

Pulse  Position  Modulation 

In  our  thesis,  we  perform  a  serial  data  communication  via  serial  port  using  RS-232  levels. 
The  data  are  transmitted  at  4800  baud  and  8  data  bits.  To  perform  serial  communication  from 
Matlab  software,  we  first  define  the  serial  port  object  and  then  configure  device  parameters  for 
communication.  Values  sent,  received,  transfer  status  and  bytes  available  are  recorded  in  Matlab 


34 


for  tracking  purpose.  We  inerease  the  buffer  size  of  Matlab  software  to  1023  bytes  (by  default  it 
is  512  byte). The  reeeived  data  is  then  loaded  to  workspaee.  We  do  the  similar  experiment  using 
Simulink  and  reeord  the  data  to  the  workspaee. 

Sinee  there  are  only  four  stationary  nodes  in  our  network,  we  use  a  statie  routing  table  that 
lists  the  routing  paths  for  all  possible  source-destination  pairs.  For  example,  we  plaee  the  four 
nodes  on  a  straight  line,  and  label  them  as  A,  B,  C  and  D  respeetively.  If  node  A  generates  a 
message  with  the  source  field  as  A  and  the  destination  field  as  D,  then  the  next-hop  is  node  B 
and  then  node  C  in  the  routing  table  for  the  (A,  D)  souree-destination  pair. 


Table  3.  Bit  loss  test  for  different  distance  of  sensors 


Distance 
between  sensors 

No.  ofbits 
transmitted 

Received  no.  of 
bits 

Bit  loss 

Bit  error  rate 

5  m 

1023 

1015 

8 

0.007 

7m 

1023 

1011 

12 

0.011 

10m 

1023 

1011 

12 

0.011 

In  Table  4,  we  transmit  1023  bits  of  binary  sequenees  data  to  the  receiver  sensor  node  at  a 
distanee  of  5  m.  This  channel  transmits  all  1023  bits;  however  it  has  bit  error  of  8  bits  in  average. 

Table  4.  Result  analysis  for  First  Hop  Txi  and  Rxi  (d=5  m) 


Bits  transmitted 

Bits  received 

BER 

I(X;Y) 

Information  loss 
H(X\Y) 

1023 

1023 

0.0078 

1.00 

0.0 

1023 

1023 

0.0087 

1.00 

0.00 

1023 

1023 

0.0078 

1.00 

0.00 

In  Table  5,  the  transmitting  data  are  exaetly  the  received  data  from  first  senor  node;  however 
we  filter  the  error  bits  that  are  ineurred  in  first  channel.  This  channel  is  also  set  for  distanee  of  5 
m  between  transmitter  and  reeeiver.  We  get  a  full  eommunieation  of  data  between  sensor  nodes 
with  bit  error  of  7  -8  bits  in  average. 
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Table  5.  Result  analysis  for  Second  Hop  Tx2  and  Rx2  (d=5  m) 


Bits  transmitted 

Bits  received 

BER 

I(X;Y) 

Information  loss 

H(X\Y) 

1015 

1015 

0.0068 

1 

0 

1014 

1014 

0.0078 

1 

0 

1015 

1015 

0.0088 

1 

0 

In  Table  6,  we  load  the  reeeived  data  from  sensor  node  3  and  transmit  exaetly  to  another 
sensor  node  in  a  relay  network.  This  reeeiver  is  again  at  a  distanee  of  5  m.  The  ehannel  has  full 
eommunieation  of  data  with  bit  error  of  approximately  7-8  bits. 


Table  6.  Result  analysis  for  third  Hop  Tx3  and  Rx3  (d=5  m) 


Bits  transmitted 

Bits  received 

BER 

I(X;Y) 

Information  loss 

mw) 

1008 

1008 

0.0079 

1 

0 

1007 

1007 

0.0079 

1 

0 

1006 

1006 

0.0089 

1 

0 

Table  7  shows  our  improvement  regarding  minimization  in  bit  loss  and  bit  error  eompare 
to  direet  transmission  between  two  senor  nodes  and  a  relay  network.  The  ehannel  distanee  is 
approximately  15  m.  The  output  result  shows  that  the  ehannel  has  bit  error  along  with  bit  loss. 
This  proves  that  transmission  loss  in  underwater  aeoustie  eommunieation  depends  on  distanee. 


Table  7.  Result  analysis  for  Txi  and  Rx4  (d=15  m) 


Bits  transmitted 

Bits  received 

BER 

I(X;Y) 

Information  loss 
H(X\Y) 

1023 

1012 

0.0185 

0.9156 

0.0844 

1023 

1015 

0.0156 

0.9419 

0.0581 

1023 

1009 

0.0215 

0.8952 

0.1048 

Table  8  shows  that  underwater  aeoustie  eommunieation  (UAC)  eapaeity  is  strongly 
dependent  on  transmission  distanee.  Using  more  than  one  sensor  as  multi-hop  network  system 
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model,  increases  the  mutual  information  and  rate  of  data  transmission  with  minimum  bit  error 
probability. 


Table  8.  Result  analysis  for  Multihop  relay  link  and  direct  link 


S.N 

Case  study 

I(X:Y) 

BER 

Information  loss  H(X\  Y) 

1 

Txi  and  Rxi 

1 

0.0081 

0 

2 

Tx2  and  Rx2 

1 

0.0078 

0 

3 

Tx3  and  RX3 

1 

0.0082 

0 

4 

Txi  and  RX4 

0.9054 

0.1853 

0.0824 
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Conclusion  and  Future  Work 

The  problem  of  an  optimum  plaeement  of  sensors  in  a  multi-hop  relay  strueture  for  long 
range  aeoustie  eommunieation  in  underwater  channel  have  been  investigated  to  maximize  the 
mutual  information  and  channel  capacity.  The  optimum  communication  frequencies  for  given 
link  distance  and  bandwidth  available  have  been  calculated  where  noise  power  is  low.  The 
deployment  of  sensors  in  multi-hop  relay  structure  for  large  transmission  range  has  been 
considered  to  minimize  bit  error  rate  and  information  loss  as  well. 

Additionally,  an  analytical  and  experimental  result  of  multi-hop  sensor  network  model  in  a 
relay  has  been  investigated  to  optimize  the  mutual  information  and  maximize  rate  of  data 
communication.  It  was  advantageous  to  accomplish  such  transmission  using  relays  since  the 
available  acoustic  bandwidth  decayed  with  increase  in  distance.  Moreover,  numerous 
applications  required  a  relay  acoustic  link  where  data  transmission  had  to  be  accomplished 
beyond  short  distances.  The  multi-hop  relay  structured  that  we  have  used  would  be  suitable  fit  to 
address  the  issues  of  power  and  energy  consumption,  and  link  delay  along  with  maximizing 
mutual  information. 

The  improvement  in  mutual  information  and  data  rate  transmission  has  been  investigated 
using  information  theoretic  tools  such  as  entropy,  conditional  entropy,  and  information  loss 
between  transmitted  and  received  data  signals.  Based  on  our  experiment  and  simulation  result, 
we  concluded  that  the  mutual  information  between  transmitted  and  received  data  signal  can  be 
maximized  by  placing  extra  sensor  nodes  in  between  a  transmission  range  which  also  decreased 
BER  and  information  loss.  The  problem  occurred  in  our  analytical  result  was  an  additional  delay 
due  to  relay  link.  However,  we  observed  that  delay  depended  on  the  rate  at  which  the 
information  is  transmitted  as  well.  This  fact  implied  an  interesting  property  of  an  acoustic  link: 
although  each  relay  introduces  an  additional  delay,  as  the  hops  become  shorter  transmission  can 
be  accomplished  faster  over  each  hop.  This  improvement  is  specific  to  the  underwater  acoustic 
environment,  where  bounds  of  mutual  information  and  bandwidth  are  dependent  on  the 
transmission  distance. 

During  our  experiment  we  had  several  issues  with  data  compatibility  with  the  software.  It 
took  several  hours  to  run  the  Matlab  programs.  Memory  buffer  size  of  sensors  was  another  issue, 
which  only  transmit  real  data  signal  through  serial  communication.  We  also  observed  that  Matlab 
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serial  communication  does  not  support  complex  signals.  We  had  problems  with  the  sensors  since 
they  transmitted  some  random  signals  by  themselves  due  to  unknown  internal  error. 

We  can  expand  our  work  for  movable  sensor  nodes  in  underwater  channel  to  address  the  issues 
with  Doppler  spread.  Along  with  a  relay  network,  we  can  address  a  multi-path  channel  model  and 
use  the  wireless  communications  knowledge  of  MRC  (Maximal  Ratio  Combination)  to  increase 
the  detection  quality  of  received  signal. 
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Appendix 

MATLAB  SIMULATION  CODE  FOR  BPSK  TRANSMISSION  OVER  AWGN  CHANNEL 


close  all; 
clear  all; 
clc; 

SNRdB=l : 1 : 12; 

SNR=10.  ^ (SNRdB/10) ; 

Bit  Length=10^6; 

BER  Simulated=zeros ( 1 , length ( SNRdB) )  ; 


%Signal  to  Noise  Ratio  in  dB 
%Signal  to  Noise  Ratio  in  Linear  Scale 
%No.  of  Bits  Transmitted 
%Simulated  Bit  Error  Rate 


%Detection  Scheme: (Soft  Detection) 

%+l  if  o/p  >=0 
%-l  if  o/p<0 

%Error  if  input  and  output  are  of  different  signs 


%%  BPSK  Transmission  over  AWGN  channel 
parfor  k=l : length ( SNR) ; 

x= (2*floor (2*rand (1, Bit  Length) ) ) -1; 
y= ( sqrt ( SNR ( k) ) *x) +randn ( 1 ,  Bit  Length) ; 

BER  Simulated (k) =length (find ( (y. *x) <0) ); %Total  number  of  bits  in 

error 

end 

BER  Simulated=BER  Simulated/Bit  Length; 

semilogy ( SNRdB,  BER  Simulated,  ' m-< ' ,  'linewidth'  ,2.0); 
hold  on 

semilogy ( SNRdB, qfunc (sqrt ( SNR) ) , ' b- ' , ' linewidth ',2.0) ; 

%Theoritical  Bit  Error  Rate 

title ('BPSK  over  AWGN  Simulation '); xlabel (' SNR  in  dB ' ) ; ylabel ( ' BER ' ) ; 
legend (' BER ( Simulated) ', ' BER (Theoritical) ') 
axis  tight 
grid 


INFORMATION  THEORETIC  TOOLS  TO  DETERMINE  MUTUAL  INFORMATION,  INFORMATION  LOSS  AND  BER 

%**********Events  Probability  Distribution*********** 

%Event  A-->  P(x=l,  y=l);  Event  B-->  P(x=0,  y=0); 

%Event  C-->  P(x=l,y=0)  ;  Event  D-->  P(x=0,  y=l); 

%Event  F-->  P(Y!  =0,  1)  No  Joint  Probability  ; (  Bit  Loss) 

A=0;  B=0;  C=0;  D=0;  F=0;  %  Events  Initialization 

n=1023;  %Number  of  transmitted  bits 

%%Adjusting  Data 
for  i=l : n; 

x(i,:)=  (-l).''i; 

end; 

sx=size (x) ; 
sy=size (y) ; 
for  i=sy:n 
y (i) =5; 
end; 
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x(x==-l)=0;  y(y==-l)=0;  %Switching (-1) --> (0) ;  for  simpler 

calculation 

z=x-y; 
for  i=l : n 

if  z(i)==0  &&  x(i)==l 
A=A+1 

elseif  z(i)==0  &&  x(i)==0 
B=B+1 

elseif  z(i)==l  &&  x(i)==l 
C=  C+1 

elseif  z(i)==  -1  &&  x(i)==0 
D=D+1 

else  F=  F+1 

end 

end 

Tot= (A+B+C+D+F) ;  %Sum  of  All  Events 

%Joint  Probabilities 
PA=  A/Tot; 

PB=  B/Tot; 

PC=  C/Tot; 

PD=  D/Tot; 

PF=  F/Tot; 

%since  binary  data; 

Pxl=sum (x (x==l ) ) /n;  Px0=sum (x (x==0 ) ) /n; 


%P (y=l, x=l) ; 
%P (y=0, x=0) ; 
%P (y=0, x=l) ; 
%P (y=l, x=0) ; 
%P(y!=l, 0) ; 


%  Estimating  the  P(y=l)  and  P(y=0); 
countl=0;  count2=0;  count3=0; 
for  i=l : n 
if  y(i)==l 

countl=  countl+1; 
elseif  y(i)==  0 

count2=  count2+l; 
else  count3=count3+l ; 
end 

end 

Pyl=countl /n;  %  P (y=l ) -->Probability  of  receiving  1; 

Py0=count2 /n;  %  P (y=0) -->Probability  of  receiving  0; 

Pyy=count3/n;  %  P (y ! =0 , 1 ) -->Probability  of  receiving  null; 


%  Conditional  Probabilities 


PylXl=  PA/Pxl; 
Py0X0=  PB/PxO; 
PyOXl=  PC/Pxl; 
PylX0=  PD/PxO; 


%P(y=l |x=l) 
%P(y=0|x=0) 
%P(y=0|x=l) 
%P(y=l |x=0) 


P(y=l,  x=l)/P(x=l); 
P(y=0,  x=0)/P(x=0); 
P(y=0,  x=l)/P(x=l); 
P(y=l,  x=0)/P(x=0); 


%  Entropy,  Joint  Entropy  and  Mutual  information  Calculation 

Hx=  (n*Pxl*log2  (1/Pxl));  %  H(x)  Entropy (transmitted  bits) 


%  H(y)  Entropy  (Received  bits); 

Hy=  (Pyl*log2 (1/Pyl) ) + (Py0*log2 (1/PyO) ) + (Pyy*log2 (1/Pyy) ) ; 


%  H(x,y)  JointProbablility 

Hxy=  (-l*PA*log2  (PA) ) - (PB*log2  (PB) ) - (PC*log2  (PC) ) - (PD*log2 (PD) ) - 
(PF*log2 (PF) ) ; 


43 


HxY=Hxy-Hy;  %Conditional  Entropy  H (x | y) =H (x, y) -H (x) ; 

HyX=Hxy-Hx;  %Conditional  Entropy  H (y | x) =H (x, y) -H (y) ; 


Ixy=Hx-HxY;  %  I(x,y)  Mutual  Information  between  Transmitter  and 

Receiver 


BIT  ERROR  RATE  SIMULATION 


%Tx=Transmitted  bits  vector 
%Rx=received  bits  vector 

%LTx=  length  of  transmitted  bits  vector 
%LRx=length  of  received  bits  vector 
Clear  all; 

%Partial  cross-correlation  from  bit  vector  Tx  and  RX : 
for  lag=l : length (Rx) -length (Tx) -1 ; 

corr (lag) =Tx. *Rx (lag: length (Tx) -1+lag) ' ; 

end; 

%correct  bit  offset. 

off=  find(corr==  max (corr) ) -1 ; 

%error  vectore  err. 

err=  Tx-Rx(off+l  :  length (Tx) +off) ; 

%Total  number  of  bit  errors. 

te=  sum (abs (err) ) ; 

%Calculate  BER. 

ber=  te . /length (Tx) ; 

PSD  CODE  FOR  NOISE 

Fs  =  80000;  %  Sampling  frequency 

load ( ' Simulink_Serial_3_2_2014NOISE .mat ' ) 

X  =  xcorr (Noise) ;  %Noise  is  the  measured  noise  data  loaded  form  workspace 
figure;plot (x) ;xlabel ( 'Autocorrelation 

Length '); ylabel (' Amplitude '); title (' Autocorrelatoin  of  Measured  Noise') 
nfft  =  2 ''nextpow2  ( length  (x)  )  ; 

Pxx  =  abs (fft (x, nfft) ). ^2/length (x) /Fs; 

%  Create  a  single-sided  spectrum 

Hpsd  =  dspdata . psd ( Pxx ( 1 : length ( Pxx) /2 ),' Fs ',  Fs )  ; 
figure;  plot (Hpsd) ; 

DESERT  STAR  SYSTEM  (VENDOR)  SERIAL  COMMUNICATION  TRANSMISSION  CODE 

clc; 

clear  all; 

% - SERIAL — COMMUNICATION — Transmission_  Code 

%  Find  a  serial  port  object. 

s  =  instrf ind ( ' Type ' ,  'serial',  'Port',  ' COM6 ' ,  'Tag',  ''); 

%  Create  the  serial  port  object  if  it  does  not  exist 
%  otherwise  use  the  object  that  was  found. 

if  isempty(s) 

s  =  serial (' COM6 ') ; 
else 

f close ( s ) ; 
s  =  s (1) 
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end 

%configure  instrument  object,  s 

set(s,  ' InputBufferSize ' , 100)  %  Serial  Port  will  hold  data  until  1024  bytes 
set ( s , ' BaudRate ',4800); 
set  ( s ,  ' DataBits  '  , 8 )  ; 

set(s,  ' FlowControl ' ,  ' none ' ) ; %  software  handshaking  is  used  to  control  data 

flow 

set  (s,  ' Timeout ',120.0) ; 

q,  g, 
o  o 

%  Connect  to  instrument  object, s 
f open ( s ) ; 

%fprintf (s,  ' ###  '  ) 

%fprintf ( s , ' R5 ' ) 

%fprintf (s, 'D' ) 

Session_l=  fread(s); 

%  Disconnect  AND  Clean  up  all  objects 
f close ( s ) ; 
delete  ( s )  ; 


DESERT  STAR  SYSTEM  (VENDOR)  SERIAL  COMMUNICATION  RECEIVER  CODE 

%  Find  a  serial  port  object. 

s  =  instrf ind ( ' Type ' ,  'serial',  'Port',  'COM4',  'Tag',  ' ' ) ; 

%  Create  the  serial  port  object  if  it  does  not  exist.. 

%  ..otherwise  use  the  object  that  was  found. 

if  isempty(s) 

s  =  serial (' COM4 ') ; 
else 

f close ( s ) ; 
s  =  s  ( 1 )  ; 

end 

0,  q, 
o  o 

%configure  instrument  object, s 

set(s, ' OutputBuf f erSize ' , 1 024 )  %  Serial  Port  will  hold  data  until  1024  bytes 
set ( s , ' BaudRate ',4800); 
set ( s , ' DataBits ' , 8 ) ; 
set ( s ,  ' Timeout ',60.0) ; 

set(s,  'FlowControl',  ' software ');  %  software  handshaking  is  used  to  control 
data  flow 

q,  q, 
o  o 

load  TREX_M10_FB5K_2004116093249_CH3 .mat 
load  trigsignal_M10  FbSk.mat 
load  TRSEQsPLIT.mat 

%  Connect  to  instrument  object, s 
f open  ( s ) 

fwrite (s,trseql2) ; 

%pause (20.0) ; 

%fwrite  (s,trseq2); 

%pause (20.0) ; 

%fwrite  (s,trseq3); 

%  Disconnect  AND  Clean  up  all  objects 
f close ( s ) ; 
delete  ( s )  ; 
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